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Abstract 

We define a class of flows with exponential kinematics termed Periodic Exponential Shear (PES) flow which involve 
periodic exponential stretching of fluid elements along with their rotation. We exhibit analytical and numerical results 
for PES flow by using the Oldroyd-B model for viscoelastic fluids. We calculate the growth in the shear and the normal 
stresses analytically as well as demonstrate that repeated application of the flow leads to stable oscillatory shear and 
normal stresses. We define a material function applicable to a periodic, unsteady shear flow and show numerically 
that this material function exhibits deformation-rate thickening behavior for viscoelastic fluids subject to PES flow. 
We demonstrate the feasibility of PES flow by presenting preliminary experimental results from a controlled-strain 
rate rheometer, using a Newtonian mineral oil. 
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1. Introduction 

Viscoelastic fluids flowing through complex geometries are often subject to intermittent flow regimes in which the 
fluid elements are both sheared as well as extended. Examples include the flow of drilling muds[ 1], which are pumped 
from the surface of an oil rig down the drill string, through the drill bit, and back to the surface via an annular gap. 
Fluid flowing through a porous rock formation has a tortuous path through the pore space, which results in periodic 
shearing and extension of the fluid. In biological systems, the flow of blood through arteries is a periodic mixed 
flow. In flow through a microfluidic cross-slot device|]2|, fluid elements aligned with the centerline are subject to pure 
elongation coupled with translation, however, fluid elements which are off the symmetry axis are subject to a rapid 
(transient) stretch and rotation. The rotating (and stretching) fluid elements eventually align with a streamline and 
are convected thereafter. The rheological response of complex fluids under combined oscillatory-extensional flow is 
poorly understood, but critically important from a processing point of view. 



Several earlier studiesll3l I4i I5L I61I7L 1 81 I9 U 1 Oil 1 1 H have focussed on so-called exponential shear of viscoelastic fluids. 
Sivashinsky et al. measured the shear stress response of polyisobutylene solutions on a parallel-plate rheometer with 
an imposed aperiodic strain-rate of the form y yx (t; A, B) = Ae Bt (A and B are variable parameters and the dot indicates 
a derivative with respect to time). Ziille et al. (0] compared the extensional and shear viscosities measured in the 
uniaxial extension and unsteady shear flow of a low-density polyethylene melt, with the extensional strain-rate kept 
equal in magnitude to the shear strain-rate. Dealy and co-workersJ5l |6[] probed the shear stress response of branched 
polymer melts subject to an aperiodic exponential strain of the form y yx (f; a, yo) = yo(<?"' - 1), with jq the strain 
amplitude and a > the growth rate of the exponential. The authors used an elastic dumbbell model with a Hookean 
spring to show that the high strains associated with such flows lead to rapid stretching of the polymer molecules 
within the fluid. Doshi and Dealy[5] concluded that exponential shear has a strong tendency to cause elongation of 
the polymer chains within the flow; however, Neergaard et al.\l\ carried out experiments with a linear polystyrene 
melt and found that the measured shear stresses in exponential shear are bounded by shear stress values measured 
during step-strain rate experiments at comparable instantaneous strain-rates. Their theoretical studies with a reptation 
model indicated that exponential shearing flow does indeed stretch entangled, linear polymer chains, but the chains are 
stretched to an equal extent during the inception of steady shear. Demarquette and Dealyfl carried out exponential 
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shear tests using polystyrene solutions on a sliding-plate rheometer equipped with both a transducer and an apparatus 
for measuring flow-induced birefringence, and measured the shear stress as well as the third-normal stress coefficient 
[^(^a.yo) = ( T -o(0 - T zz(0)/1 / ?,i( f )] as a function of the time. Their measurements showed that the shear stress in 
exponential shear fell below the linear viscoelastic prediction. 

Venerus|9] utilized an aperiodic strain of the form y yx (t; a) — 2 sinh(af), and measured the first normal stress differ- 
ence [Ni(t; a) = r xx {t) - T vv (f)] on a shear rheometer with cone and plate geometry, using a low-density polyethylene 
melt. He compared experimental data from constant strain-rate planar elongational and exponential shear flows, 
and concluded that exponential shear is incapable of generating stresses comparable to planar elongational flows, 
due to the non-zero vorticity in the flow (see Section 2 below). Kwan et al. [10] carried out Brownian dynamics 
simulations of freely jointed bead-rod and bead-spring chain models subject to an aperiodic strain-rate of the form 
y yx (t; a) — 2a cosh(af), and compared their results with simulations of planar and uniaxial extensional flows. Graham 
et a/. [Ill] simulated the "pom-pom" molecular model for branched polymer melts under aperiodic strain-rate of the 
form y yx (t;a) — ae at and com pare d their results with the experimental data due to Ziille et al.[4] and Venerus|0]. 
Several authors H 03 QJ 1411 have defined and measured appropriate material functions for exponential shear, 



which we discuss below in Section 3.4. 

In this paper, we introduce a non-linear rheometric deformation that we term Periodic Exponential Shear (PES) 
flow, which may be used for the investigation of material response to flow fields in which rapid stretching and rota- 
tion occur simultaneously and periodically. We evaluate analytical and numerical results using the Oldroyd-B model 
for viscoelastic fluids and define an appropriate material function for PES flow, which shear-thickens (i.e., the ma- 
terial function for the flow grows faster than the linear viscoelastic prediction, as a function of the shear-rate). We 
demonstrate the utility of such a flow by showing preliminary experimental results from a shear rheometer, using a 
Newtonian mineral oil. 

2. Definition 

A general flow kinematics which incorporates periodic exponential stretching and rotation of the fluid elements, 
may be defined as: 

Jyxit) = 7o sinh[a(f - f )] for f < t < h, 

Yyx(t) = To sinh[a(fi - to)] - J\ sinh[fo(f - fi)] for t\ < t < t\ + \ sinh -1 [ — sinh[a(fi - f )] |, (1) 

b \7\ ) 

with b + 0, j\ ^0. Here a and b are the growth rates (units s ) of the exponentials, with strain amplitudes yo and 
j\. The half-periods of this flow [fi - fo, \ sinfT 1 sinh[a(fi - fo)])] are unequal in general. For simplicity here, we 
choose identical growth rates a = b and strain amplitudes yo = Ti in Eq. (Q}, and define a periodic reversing (or PES) 
flow as: 

Jyxit) = yo sinh[a(f - f )] for t <t<t u 

7yx(t) = ro sinh[a(fi - f )] - 70 sinh[a(f - fi)] for t x <t <2t { - t . (2) 

The half -periods of this derived flow are equal [t\ - to, t\ - fo] and the time period for one cycle of deformation is T — 
2(ti -fo). The parameters of PES flow are the exponential growth rate a, the strain amplitude 70 as well as the initial (fo) 
and half -period time (fi). In Fig. QI a X we plot the normalized strain as a function of the time for varying exponential 
growth rates a over one cycle of oscillation with 70 = 1 , fo = 0s, f 1 = 1 s. In Fig. [TJbX the periodic nature of PES flow 
is demonstrated for a = 3.1i _1 ,7o = 1. The velocity field in an unsteady shear flow is given by v = (y yx (t)y, 0, 0). PES 
flow is rotational, and the vorticity [a> = V x v = -y yx (t)i] grows exponentially with time and oscillates. By contrast, 
planar elongational flow [v = (ax, -ay, 0)] and uniaxial extensional flow [v = (ax, -ay 12, -az/2)] are irrotational, and 
fluid elements within such flows are subject to pure extension alone. We note that a frequently used flow-classification 
scheme (which does not take into account the vorticity of the flow) based on the invariants of the Finger strain tensor 
C~'(f, fOSH] yields the obvious result that PES may be classed as an unsteady shear flow and therefore the first 
and the second invariants of the tensor are equal (Ii = tr[C~ l (t, t')] - h = fr[C(f, f')], here tr indicates the trace of 
the tensor), as for the case of planar elongational flow. The flow kinematics defined above ensure continuity of the 
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strain over each cycle, but the strain-rate (and therefore the vorticity) is discontinuous in both sign and magnitude at 
the half -period time t — t\. Our flow kinematics is experimentally realizable on a shear rheometer (see Section 3.5 
below), and directly comparable with the more traditional exponential shear flow in the first half -period. 




Figure 1: (a) PES flow over one cycle of oscillation with yo = 1, to — O.v, fi = Is and varying exponential growth rates a: circle (a = Is '), star 
(a = 2s~'), square (a = 4.v~'). (b) Periodic nature of PES flow demonstrated for a = 3.1s ,70 = 1. 



3. Results 



For our analytical and numerical studies, we use the well-known Oldroyd-B model[15], which serves as a useful 
representation for the flow of dilute polymeric liquids: 



t + Ait w = T] (y m + A 2 y (2) ). 



(3) 



Here t is the stress tensor, y is the strain tensor, 770 is the zero-shear-rate viscosity, A\ and A 2 are representative 
time-scales for the polymer and the solvent, respectively. The subscripts indicate the order of the upper-convected 
derivative, with cT(i) = ^ + v • Ver - {(Vv) r ■ cr + cr ■ (Vv)} for arbitrary operator cr. The Oldroyd-B equations for an 
unsteady shearing flow are given by: 



(4) 



1 + M j t I r xx (t) = 2A,y yx (f)T yx {t) + 2 m A 2 yl x {t) (5) 

with T yy = t zz = 0. The shear and the normal stress are assumed to be specified at some instant of time t = fo, 
r yx (to) = t^ x and T x .t(fo) = t° xk . As discussed below, the combination Q = A\a arises frequently in the solution to the 
above system of equations and we regard it as an effective dimensionless growth rate for PES flow. The numerical 
value of Q does not determine the flow dynamics uniquely, since widely differing values of A\ and a lead to different 
shear and normal stresses, but may still represent the same numerical value of Q. We note that the Oldroyd-B equations 
are non-singular for Q — {0.5, 1}, but under PES flow, the shear stress T yx (t) is found to be singular for Q — 1 and 
the normal stress r xx (t) is found to be singular for Q = {0.5, 1}, as shown below. Hence, we solve for the different 
cases separately, in what follows. For our numerical work, we simulate Eqs. © and (|5]l using Matlab Ver. 7.9 (The 
Mathworks Inc., United States). 
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3.1. Shear stress T yx {i) 

The solution of the linear first-order ordinary differential equation [Eq. (HJl] for to < t < t\, Q + 1 with the 
strain-rates calculated from the PES flow strain profile [Eq. (0] is: 



T yx {t) 



o T]oyoa(l - a z AiA 2 ) 

Tyx + 



(A ia ) 2 - 1 



(^ia) 2 - 1 
T]oyQa 2 A 2 
'(Aia) 2 - 1 



(6) 



\A\a cosh(a(f - Jo)) - sinh(a(f - ?o))}- 



The solution for t\ < t < 2t\ — to, Q + 1 is: 
77 y fl(l -a 2 AiA 2 ) 



Ty X (t) 



(A ia ) 2 - 1 



,-(«-«i)Mi 



(A ia ) 2 - 1 
T] y a 2 A 2 



{Aia sinh(a(f - fi)) - cosh(a(f - f;))} 



(7) 



(A ia ) 2 - 1 



\A\a cosh(a(f - t\ )) - sinh(a(f - f i ))}. 



Here r' equals the value of T yx (t\) obtained from Eq. (0, which ensures continuity of T yx (f) at t = t\. In Fig. |2|a), 
Eqs. (O and (O are plotted as red curves, and shown to agree with the numerical solution (solid blue circles) of the 
Oldroyd-B equations over one cycle of oscillation for the representative parameter values a = 3s~ l ,yo = 1,770 = 
1 Pa.s., Ai = 2s, A2 = Is. In Fig. |2|b), we plot the shear stress as a function of the normalized time for 10 cycles 
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Figure 2: (a) A plot of the analytical solution (red curve) for the shear stress r yx (t) as a function of the normalized time tjX\ and a comparison with 
the numerical solution (solid blue circles) of the Oldroyd-B model for a = 3s~' ,70 = 1, t]q = 1 Pa.s., Ai = 2s, A2 = Is (here Q = 6). (b) A plot of 
the shear stress T yx (f) as a function of the normalized time t/Ai with the same parameters as above, for 10 cycles of oscillation. 



of oscillation with the same set of parameter values. For Q = 1, Eq. © is solved with a = l/A\, A\ + using the 
strain-rates calculated from Eq. (f2)- The solution for to < t < t\ is: 



_ 77oyoa(l + aA 2 ) 

• V Y a 



T yx (t) = 

The solution for t \ <t<2t\— to is: 

77oroa(i +C1A2) 



e -^ ) + Vo7 a\l-aA2) (t _ to)e _ a(t _ to) 



j]oy a(l + aA 2 ) a(t _ t ) 

H : e . 



Ty x (f) = 



T' + 



4 



-a(t-h) 



rjojoa 2 {\ -aA 2 ) 



(t-h)e 



-a(t- h ) _ W0«(1 +^2) „(,-,,) 
4 



(8) 



(9) 



with Ty X equal to the value of T yx (t\) obtained from Eq. © to ensure continuity of the shear stress at t = t\. We 
note that the Newtonian limit A\ — A 2 — > [a — » 00] is not defined for Q = 1, since the above expressions diverge 
in the limit. It appears that the singular case Q — 1 may be unphysical. In Fig. [3ja), we plot Eqs. (O and (O for 
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Figure 3: (a) A plot of the analytical solution (red curve) for the shear stress r vr (f) as a function of the normalized time tjX\ and a comparison with 
the numerical solution (solid blue circles) of the Oldroyd-B model for a = 0.5s _I , yo = 1, f/o = 1 Pels., Ai = 2s, Ai = Is (here Q = 1). (b) A plot 
of the shear stress as a function of the normalized time with the same parameters as above, for 10 cycles of oscillation. 



the representative parameter values a = 0.5s , yo = 1, 770 = 1 Pa.s., A\ = 2s, A2 = Is as a red curve, and compare 
with the numerical solution (solid blue circles). The shear stress is plotted as a function of the normalized time for 10 
cycles of oscillation in Fig. [3jb). 

The shear stress, after an initial transient, approaches a stable oscillatory state (which has been termed alternance $X2\ ) 
in both Figs. |2jb) and [3b) and we now find the analytical expression for the peak shear stress at an arbitrary time 
f = n (n is a positive integer). For the case @ + 1, we plot the shear stress as a function of the time for 10 cycles 
of oscillation in Fig. Ufa) and refer the reader to the figure for the steps that follow : In the first cycle of oscillation 
[fo = 0s, t\ - Is], the numerical value of the shear stress at t — Is is T yx (l) = (r° + a)e~ 1 ^ 1 + (J3/-io + yn\), with 

a = ' 1 ° 7 °^^_ A 1 1 ' i2 \ /i = a[j!$~[ >y = (tj'ap-i ' A*o = ha sinh(a) - cosh(a), p,\ = /liacosh(a) - sinh(a). The numerical 
value of the shear stress at t = 2s is r yx (2) = (t* - a)e~ 1/il - (flfio + y/Ji) = (t° + a)e~ 2lli + (fifxo + yn\)e~ llM - 
oier^lM _ (jg^y + yj_i x ) 7 on substituting for tL = T yx (l). The procedure may now be repeated iteratively for t — 3s to 
obtain r„(3) = (t° + a)e~ 3/Ai + (J3/j. + y/Ji)e' 2/A ' - ae~ 2/A] - (J3fi + y/u{)e~ 1/M + ae~ l/Ai + (J3fi + ym). By inspection 
of the foregoing expressions, we write below the expression for the numerical value of the peak shear stress at t — n 
(n is a positive integer): 

n 

T y M = r>- nMl + [a + e m > (J3m> + Yfii)] Y J (-l) m+n e- m/A > . (10) 

m=l 

From Eq. ( TTOt we find that T yx (l9) - 0.4382 Pa for the representative parameter values a = li -1 , yo = 1, J70 = 1 Pa.s., 
A\ = 2s, A2 = Is. The calculated value is found to agree with the numerical value of r yx (l9) in Fig. [4] In Fig. 0b), 
we plot the shear stress as a function of the shear strain for the same values as stated above, and note the approach to 
a stable oscillatory state, with a nearly elastic response at long times. For Q — 1, using a similar iterative procedure, 
we find that the numerical value of the shear stress at t — n (n is a positive integer) is given by: 

n 

Tyx (n) = (r° yx - a)e- an + ae'"^ + a(-l)"(l - e a ) + a/3 ^(-l)'" + " e - flm , (11) 

m= 1 

with a — riayaa( ^ +aAl ^ p — woflQ-a-^) 
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(19,0.4382) 




(here Q = 2). The coordinates of the peak shear stress for the tenth cycle of oscillation are indicated, for purposes of comparison with the 
analytically calculated value (see text for discussion), (b) A plot of the shear stress T yx (t) as a function of the shear strain y vv (f) for the same 
parameter values as in (a). 



3.2. Normal stress T xx (t) 

The normal stress may be calculated by substituting Eqs. © and <j7j> into Eq. <(5j and solving the resultant linear 
first-order ordinary differential equation. For to < t < t\, Q + {0.5, 1} the normal stress is found to be: 



?7o(yofl) 

"4(^ia) 2 -l 



A 2 + 



r xx (t) = C ie -' /A > + 2 7o e- ( '-"> )/A ' sinh(a(f - f )) 
A l (a 2 A l A 2 - 1) 



r° yx + 



J] -y a(l -a 2 A\A 2 ) 



(Aia) 2 - 1 



(12) 



(Aid) 2 - 1 



[4(/lio) z - 1 + 2A x a sinh(2a(f - f )) - cosh(2a(f - f ))] 



7]q(Ai - A 2 )(Axd)(jod) 
'[(A ia ) 2 -l][(4(A ia ) 2 -l] 



[2A\a cosh(2a(f - fo)) - sinh(2a(f - to))], 



with Ci = e' o//i 
normal stress is found to be 



, no(7oaf(2-4(Aia) 2 ) 



I, A i (a 2 A l A 2 -\) \ 2 ?0 (^ 2 -^i)(i 1 a) 2 ( yo a) 2 
\ A1 ^ (Ma) 2 -l )^ [(M") 2 -\][4{A,a) 2 -\} 



m(yo a ) 2 

~4(A ia ) 2 - 



A 2 + 



T xx (t) = C 2 e-' / *> -2y e-^ A ' 
A l {a 2 XiA 2 -l) 



sinh(a(f - fi)) 



For fi < t < 2t\ - to, Q + {0.5, 1} the 
! j]ojoa(\ - a 2 A X A 2 ) 



(Aid) 2 - 1 



(13) 



(Aid) 2 - 1 

j]o(A 1 - A 2 )(Aid)(yoa) 2 



[4(Aid) 1 - 1 + 2A\a sinh(2a(f - fi)) - cosh(2a(f - fi))] 



[(A ia ) 2 - l][(4(A ia ) 2 - 1] 



[2A\a cosh(2a(f - 1\)) — sinh(2a(f - t\))], 



with C 2 — e 



ti/M 



n (y a) 2 (2-MAia) 2 ) 



(a 2 



A l (a I A l A 2 -l) 



) 



lUM-AQiA^fiyoa) 2 
[y,fl) 2 -l][4(iia) 2 -l] 



. Here t\ x is the value of T xx (t\) obtained 



4(,lia) 2 -l V 12 T (^ia) 2 -l 

from Eq. (TTZb . which ensures continuity of the normal stress at t = t\. As for the case of the shear stress, special 
solutions for the normal stress may be found for the singular cases Q - {0.5, 1}, but are not discussed here. In Fig. 
|3a), we plot Eqs. (fT2l and (fT3l l for the representative parameter values a = 3s~ l , jo = 1, r\o — 1 Pa.s., A\ = 2s, 
A 2 = 1 s as a red curve and compare with the numerical solution (solid blue circles) of the Oldroyd-B model. The 
evolution in the normal stress with the same parameter values is plotted as a function of the normalized time for 10 
cycles of oscillation in Fig. [2b). The normal stress is shown to approach a stable oscillatory state after a complex 
initial transient response and the analytical expression for the numerical value of the peak normal stress may also be 
found via an iterative procedure, as outlined above for the shear stress. 
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Figure 5: (a) A plot of the analytical solution (red curve) for the normal stress r xx (t) as a function of the normalized time t/A\ along with the 
numerical solution (solid blue circles) of the Oldroyd-B model for a = 3s _I , yo = 1, 770 = 1 Pels., Ai = 2s, A% = Is (here Q = 6). (b) A plot of the 
normal stress T xx (i) as a function of the normalized time t/Ai with the same parameters as above, for 10 cycles of oscillation. 



We note that for the Newtonian limit of the Oldroyd-B model: A\ = X% — > 0, Eqs. (|6j, (0 reduce to r yx (t) = 
77070a cosh(q(f — fo)) = rioy yx (t), while Eqs. (fl2t . ([T3T i reduce to t xx = 0, which are consistent with the anticipated 
expressions [121 for a Newtonian fluid. 

3.3. Derived quantities 

In oscillatory shear, the kinetic energy dissipated per unit volume per cycle of imposed strain oscillation is given by 
e = j T yr dv y Al5\\ . For PES flows, we use Eqs. (O and (|7]i to calculate the dissipation rate for Q + 1, to < t < 2t\ — fo: 



e = 



(Aia) 2 - 1 



is 277oyoa(l - a 2 AiA 2 ) 

i T yx ~ T yx> + 



(A ia ) 2 - 1 



-(f!-f )M 



(A\a sinh(a(fi - f )) + cosh(o(fi - f ))) - l] (14) 



2[{Aia) z - 1] 2[{AiaY - 1] 

In the Newtonian limit A\ — A2 — > 0, the above expression reduces to e — ^j~[2a(ti - to) + sinh(2a(?i - fo))]- 

A measure of the average orientation of polymer chains within a viscoelastic fluid flow is given by the extinction 
angle ^(f) = ^tan~ l [2T yx {f)IN\{t)\ |7]. In Fig. |6| a), we plot x(t) as a function of the normalized time for the parameter 
values a = 0.001s -1 , yo = 1, ^0 = 1 Pa.s., A\ = 2s, A 2 = Is. The limit a — > is found to be equivalent to the 
Newtonian limit of the Oldroyd-B model, and the graph of the extinction angle alternates between x ~ an d 
X ~ —it I A every half-cycle. In Fig. |6fb), we plot the extinction angle as a function of the normalized time for similar 
parameter values as chosen for Figs. |2}t>) an d EJb). We find that after an initial transient, the polymer chains are 
stretched and rotated periodically about the flow direction as a result of the vorticity in the flow. 

3.4. Material function for PES flow 

In earlier studiesflEi 12, 14] several different material functions have been proposed for exponential shear. Ziille 



et al. |01 used a material function r\z (where the subscript Z indicates "Ziille") which varies with the choice of the 
coordinate system employed and does not account for contributions from the normal stresses: 

T (0 

T] Z (t) = (15) 

a 

Although r/z shows shear-thickening behavior for the exponential shear of viscoelastic fluids, Samurkas et al.^ 
pointed out that the material function also shear-thickens for a Newtonian (purely viscous) fluid, on account of the 
exponential growth of the shear stress, so unbounded shear stress growth does not imply shear-thickening behavior of 
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Figure 6: (a) A plot of the extinction angle ^(/) = 
PES flow with the parameter values a = 0.001.s~ 



(0.5)tan 1 [2r >v (f)/r„(f)] as a function of the normalized time t/Ai for 10 cycles of oscillation of 



,T0 : 



l,i/o = 1 Pa.s.,/li 



■ 2s, /l? 



l.v (here Q = 0.002). (b) A plot of the extinction angle as a 



function of the normalized time for 10 cycles of oscillation of PES flow with the parameter values a - 
(here Q = 6). 



3s ,yo = l,i/o = 1 Pa.s.,/ii = 2s, A2 = Is 



the fluid just the exponential increase of the deformation rate. For steady shear flows, Jobling and Roberts[13(|, and 
independently, Lodge and Meissner lll4ll proposed a shear viscosity tjlm (subscript LM indicates "Lodge-Meissner") 
which is independent of the choice of any special coordinate system used to define the flow, and is given by the ratio 
of the difference of the first two components of the principal stress tensor divided by the difference of the first two 
components of the principal strain-rate tensor, which may be rewritten in terms of the first-normal stress difference 
and the shear stress: 



2\y yx 



(16) 



Doshi and Dealy|5] suggested that Lodge and Meissner's definition of the shear viscosity be applied to the case of 
unsteady shear flow, with the difference of the (time-dependent) principal strain-rates in the denominator: 



riDD(t) = 



^(0+4x2,(0 



(17) 



2\j yx {t)\ 

where the subscript DD indicates "Doshi-Dealy". 

We note that Eqs. (1151 1. ( TTol l. and (fTTb are clearly inappropriate for the case of PES flow, which is a periodic 
unsteady shear flow with a discontinuous strain-rate at the half -period time t \ . We propose a frame-invariant material 
function, based on the Lodge-Meissner viscosity, which we call the Generalized Lodge-Meissner viscosity t]glm, 
suitable for any periodic unsteady shear flow, calculated over an oscillation cycle: 



1]GLM = 



§ ^N\{t)+Ar 2 yx (t) dt 
2§\y yx {t)\ dt 



(18) 



The Generalized Lodge-Meissner viscosity is calculated by integrating the difference of the first two components of 
the principal stress tensor over an oscillation cycle and dividing by the integrated difference of the first two components 
of the principal strain-rate tensor over the same cycle. We note that for a Newtonian fluid (defined by T yx {f) = J]j y _ x (t) 
and Ni(t) = 0), the Generalized Lodge-Meissner viscosity reduces to the anticipated result t\glm = J]- For PES 
flow, upon substituting for the strain-rate in the denominator of Eq. (II St . we obtain the PES viscosity rjp (where the 
subscript P indicates "PES flow") calculated over a single oscillation cycle: 



4yo sinh[a(?i 
8 



to)] 



(19) 



We plot Eq. ( fT9l with the representative parameter values yo = 1, 770 = 1 Pa.s.,/li = 2s, A2 = Is and varying 
exponential growth rates a in Fig. 0a). For convenience, we plot the calculated value of rjp over each cycle, at the 
mid-point of the cycle. We note that the PES viscosity always approaches a steady-state value and is continuous. In 
Fig. |7jb), we plot the calculated steady-state value (after 100 cycles of oscillation) of the normalized PES viscosity 
rjp^sMo as a function of the normalized exponential growth rate aT (here T = 2{t\ - to) is the time-period) for 
yo = 1,770 = 1 P a - S - an d varying relaxation times A\,A2 = Aj/2. We note that the steady-state PES viscosity grows 
faster than the linear-viscoelastic prediction (which corresponds to rjpss /% = 1 ), indicating shear-thickening behavior 
of the viscoelastic fluid. By contrast, Neergaard et concluded from numerical studies with a reptation model that 
(aperiodic) exponential shear does not exhibit shear-thickening "in any practicable experiment for linear, entangled 
polymer chains in shear flow". 




Figure 7: (a) A plot of the normalized PES viscosity )//>/?/o as a function of the normalized time t/A\ for 20 cycles of oscillation with yo = 1, tjq = 
1 Pa.s,/li = 2s, Xj = Is and varying exponential growth rate a with the legend: circle (a = 0.5.v~'), square (a = Is -1 ), diamond (a = 1.5s _1 ), 
star (a = 2.v~'). (b) A plot of the steady-state normalized PES viscosity rjpss /vo as a function of the normalized exponential growth rate aT (here 
T = 2(t[ - tg) is the time-period of one flow-cycle) for varying relaxation times Ai, A2 = A[ /2 with the legend: star (A[ = Is), circle (Ai = 2s), 
square (Ai = 4s) and yo = l,r/o = 1 Pa.s. 



3.5. Experimental results 

Our experiments were carried out at room temperature on a strain-controlled Advanced Rheometric Expansion 
System rheometer (ARES-LS8 18409, TA Instruments, United States). We used a cone-plate assembly with a cone di- 
ameter of 50mm (cone angle=0.04rad). The rheometer permits direct acquisition of dc voltage signals from the torque 
transducer, the optical encoder, and the transducer normal force sensor through BNC connectors in the rear panel of 
the instrument. These unprocessed voltage signals are not noise-filtered or corrected for inertia and compliance of the 
torque transducer. Data was acquired at 16-bit resolution through an analog-to-digital card (NI PCI-6032E, National 
Instruments, United States) coupled with a Labview (National Instruments, United States) code at a sampling rate of 
10 3 points per second. The acquired signal was filtered for noise using a Savitzky-Golay filter and calibrated to find 
formulae which were used to convert the voltage values to quantities of physical interest. The calibration curves used 
were y = 0.035x (x in volts, y in Newton-metre) for the torque, y — Q.lx (x in volts, y in radians) for the deflection 
angle and y = -4.05jc (x in volts, y in Newtons) for the normal force. The values of the stress (in Pascals), the strain 
and the first-normal stress difference (in Pascals) were calculated from the torque, the deflection angle, and the normal 
force respectively, using conversion factors appropriate to the measuring system geometry and torque transducer em- 
ployed. The rheometer can be easily programmed to accomplish exponential shearing flow by supplying the desired 
strain profile to the proprietary software used to control the instrument. However, due to a hardware limitation of 
the rheometer, continuous oscillation cycles of PES flow could not be achieved: The rheometer accomplishes two 
complete oscillation cycles of PES flow, after which there is a delay of approximately 25 s (during which the specified 
waveform is sent to the instrument's memory) before the next two cycles commence; which is why we show only two 
continuous cycles in the plot that follows. 
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For our tests we used a viscous Newtonian mineral oil N1000 (Cannon Instrument Company, United States) with 
the stated zero-shear-rate viscosity 770 = 2.01 Pa.s at 25°C. In Fig. [3a), we plot two cycles of oscillation of the 
imposed strain as a function of the time for jq — 1, a — 3.2.S -1 . The measured shear stress response for N1000 oil is 
plotted as a function of the shear strain (red curve) in Fig. 13b). The Lissajous plot shows that the rheometer has a lag 
time of approximately 0.05s at startup of motion of the plate. For purposes of comparison, the shear stress calculated 
numerically from the Oldroyd-B model with parameter values a = 3.22s _1 , 770 = 1.48 Pa.s., jq = \,A\ = X% = 10~ 4 
(which approximate the Newtonian limit of the model) is also shown (blue curve). The measured normal force for the 
Newtonian N1000 oil was found to be a small noisy signal, and is not shown here. 




Figure 8: (a) A plot of the imposed shear strain y yx (f) as a function of the time t for two cycles of oscillation of PES flow with a = 3.2s ,yo = 1. 
(b) A comparison of the response shear stress r,, v (f) for two cycles of oscillation as a function of the shear strain y yx (t) for the Newtonian mineral 
oil N1000 (red curve) with the numerically calculated shear stress (blue curve) having model parameter values: a = 3.22s~',7/o = 1.48 Pa.s., 
yo = 1,^1 = /I2 = 10~ 4 , which approximate the Newtonian limit. 



4. Conclusion 

To summarize, we have introduced a new rheological test protocol consisting of a periodic, exponential shear 
flow, and discussed the response of a viscoelastic fluid such as a dilute polymer solution that is well-described by the 
Oldroyd-B constitutive equation. Exact analytic expressions for the shear and the normal stress have been calculated 
and compared with numerical computation of the model. An appropriate material function for this periodic, oscillating 
flow has been defined and shown to exhibit deformation-rate hardening. 

It is known|l]] that fluid invasion into the well-bore of an oil rig weakens the rock formation, and can cause a collapse. 
Improved well-bore strengthening may be achieved by the use of drilling fluids which shear-thicken while flowing 
through the drill bit. The observation (or lack thereof) of shear-thickening in the drilling fluid could potentially be 
related to the control of fluid invasion near the well-bore. The PES flow kinematics defined in this paper may serve 
as a basis for distinguishing between different formulations of drilling fluids, by studying their rheological response 
to a well-characterized large-strain oscillating flow. We hope that our comprehensive study of a periodic flow which 
incorporates both rotation and rapid stretching of the fluid elements is found to be of utility within laboratory and 
industrial settings. 
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